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1 Introduction 



In classical mechanics a lot of integrable systems on so(4) and e(3) algebras are related with 
the integrable spin chains. For instance pQ, different Gaudin magnets coincide with the Euler, 
Lagrange, Neumann and Clebsh systems on e(3) and with the Manakov, Steklov systems on 
so(4) . The separated variables for all these models may be derived from the separation of 
variables for the XYZ Gaudin model 

Some degenerate cases of the Heisenberg spin chain are connected with the Goryachev- 
Chaplygin top auxiliary symmetric Neumann system and Kowalevski-Goryachev-Chaplygin 
top |3|. The separated variables for all these models may be derived from the separation of 
variables for the XYZ spin chain [Hj. 

In this note we consider XXX Heisenberg magnet with boundaries [BJ [7] using the Lax 
matrix for the standard two-site XXX Heisenberg magnet 



T(A) 



A - s 3 + iSi s 1 + is 2 
si — is 2 A + s 3 + i8i 



A - t 3 + i5 2 t\ + it 2 
ti - it 2 A + t 3 + id 2 



XI) 



Here Si are numerical shifts of the spectral parameter A. Dynamical variables Sj,£, are coordi- 
nates on so(4) = so(3) © so(3) with the following Lie-Poisson brackets 



where e^k is the totally skew-symmetric tensor. 

Matrix T(\) p.lj) defines representation of the Sklyanin algebra 



1.2) 



[r(A-i/),T(A)T(i/) 



1.3) 



on generic symplectic leaves of so(4). Here we use the standard notations T(X) = T(A) <g) 
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/ , T(v) = I (g) T{y) and r-matrix has the form 



r(A — v) 
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:i-4) 



Applying the standard machinery jHl U\ to T(A) one gets another integrable systems with 
cubic and quartic additional integrals of motion. As a matter of fact, the construction leads 
automatically to the separated variables. 



2 Cubic integrals of motion 

The main property of the Sklyanin algebra (jl.3|) is that for any numerical matrix K, coefficients 
of the trace of matrix /CT(A) give rise the commutative subalgebra 

{tr/CT(A),trO»} = 0. 
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All the generators of this subalgebra are linear polynomials on coefficients of entries Ty(A), 
which are interpreted as integrals of motion for integrable system associated with matrix T(X). 
For instance, representation (jl.ljl generates one linear and one quadratic integrals of motion 
in variables Sj,ij and the corresponding integrable system is equivalent to a special case of 
Poincare system jSJ. 

According to |Zj , we can construct commutative subalgebra generated by quadratic polyno- 
mials on coefficients of I^(A). Let us introduce the matrix 

f(X) = JC d (X) T(A) , (2.5) 

where T(X) is given by and 

^w = (m + 1, M o +ao )- < 2 ' 6 ' 

Here a,k,bk are arbitrary numerical parameters and Aq depends on dynamical variables 

A = a 1 (is 2 + it 2 -h- si) - 6i(si +t x + is 2 + it 2 ) - s 3 - t 3 . 

We can say that dynamical matrix K d ()2.6|) describes some dynamical boundary conditions. 
The trace of T(X) 

f(A) =trf(A) = A 3 + hX + I 2 

give rise the commutative subalgebra of the Sklyanin brackets (II. Hj) . Quadratic and cubic in 
variables Si,ti polynomials I\ y2 are integrals of motion for some integrable system on so(4). 

To compare this system with the known examples of integrable systems on so(4) we intro- 
duce two vectors J = (Ji, J 2 , J3) and y = (2/1,2/2, 2/3) with entries 

Hi = x(si-ti), Ji = Si + U, (2.7) 

which satisfy to the following Lie-Poisson brackets 

{ Jj , Jj } = Sijk Jk 7 {^i ! Uj } = ^rjfc 2/fc , {2/1 1 2/ j } = >^^ijkJk ■ (2.8) 

Because physical quantities 2/fc, ^fc should be real, parameter x 2 must be real too and alge- 
bra ()2.8|) is reduced to its two real forms so(4, R) or so(3, 1,R) for positive and negative x 2 
respectfully. The corresponding Casimir elements are equal to 

C„=x 2 |J| 2 + |y| 2 , C=(y,J). (2.9) 

Let (y, J) and y x J stand for the inner vector product and for the vector cross product 
respectively. In variables ()2.7j) the Hamilton function is equal to I\ up to constants 

H = 2/ 1 + ^ + 25 1 5 2 = (2.10) 
2x z 

= I J| 2 + (a, J)(b, J) + x^b, y x J) - (b, (5 1 - ^x^y + (5 1 + 5 2 ) j) + 2(c, J), 
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where numerical vectors are 

a=^0, 0,2iJ, b = (i(ai + bi), a\ — &i, zj, c = (a + b , -i(a - b ), . 
Additional integrals of motion K looks like 
K = -4zJ 2 = (b,J)[2|J| 2 + x- 1 (a,y x J-(5 1 -5 2 )y + x(5i + <5 2 )j) -x- 2 ^-45i5 2 

- 2x-\5 l - 5 2 )(c, y) + 2x- 1 (c,y x J) + 2(5 1 + 5 2 )(c, J) (2.11) 

Integrals of motion H and K are defined up to canonical transformations. 

Suppose that Hamiltonian function has to depend on the third component of the vector 
y x J only. The Hamiltonian ()2.10|) has such form after rotation y — > Uy and J — > UJ on the 
following Euler angles 

7r In di — In b] , it , . / — ~ . % , . 

0=--2 L_ ^, = 0, = --zln(-^ + 2 v / aA) --ln(4oi6i + l). 

This rotation acts on the numerical vectors a, b and c in the following way 

a = [/-i a = (o, 2 V / T^c 2 ", 2zc) , b = U' l h = (o, 0, -wT 1 ), 



E=IT'c=i (eft -^^/J 

where 

1 2«(ai6 - a 6i) fl -2(ai&o + a 6i' 

c = =, a = — , p - 



V4ai6i + 1 ' v /a 1 fe 1 (4a 1 6 1 + 1)' v /a 1 fo 1 (4a 1 b 1 + 1) ' 

Substituting a, b and c instead vectors a, b and c into (|2.1U|) and ()2.11|) one gets integrals of 
motion after rotation. The Hamilton function H (j2.1(J|) after rotation and renormalization 

H 



H = = c{Jl + J 2 - Jl) - 2V1 - c 2 J 2 J 3 + — (y 2 Ji - V1J2) + ctJx + ii-h + jJ 3 + Sy 3 . 

V4ai6i + 1 ix 

(2.12) 

depends on five essential parameters c,a,/3 and 

7 = -1(61 + 5 2 ) +4— — — — , 5 = — . 

4axOi + 1 ix 

It is a real function on so(3, 1,M) with negative x 2 only. 

According to jHUZj, the separated coordinates gi )2 for (|2.12|) are zeros of the polynomial 

Tn (A) = (A-gi)(A-g 2 ) = 0, (2.13) 
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whereas the conjugated momenta are equal to 

p k = -i In T 2i (q k ) - ln(ai q k + <k>) ■ (2.14) 

We can prove that q k ,p k are Darboux variables using (|2.13H2.l3)l and brackets (|1.3jl . 
By definition the generating function of integrals of motion is 

f (A) = trace T(A) = A 3 + ^A - 7 2 = 

= (A + A ) Tn (A) + (aiA + a ) T 21 (A) + (61A + b ) T 12 (A) 
Substituting X = q k into this equation one gets two separated equations 

9fe + h qk + h = exp(zpfc) + det T(g fe ) exp(-zp fc ), k = 1, 2. 
Here we took into account that Tn(q k ) = and Ti 2 (q k ) = detT(q k ) T^iqk). 



3 Quartic integrals of motion 

According to [U], we can construct another commutative subalgebra generated by quadratic 
polynomials on coefficients of Tjj(A), which are integrals of motion for another integrable system 
associated with the same matrix T(A) (jl.lj) . 

Let /C±(A) be generic numerical solutions of the reflection equation 

^ _ ( b 3 X + a (b 1 +ib 2 )X\ _ / a 3 A + /3 (01 + ia 2 )A \ . . 

+ I (61 — «6 2 )A -b 3 A + a y ' I (ai — ia 2 )A -a 3 A + /3 y , K 1 



The Lax matrix for the two-site XXX Heisenberg magnet with boundaries has the form 

L(X) = /C + (A)T(A)/C_(A) (\ I ) T*(-A) ( ^ J 
Trace of 7(A) 

r(A) = trL(A) = -2(a,b)A 6 - I X A 4 - / 2 A 2 - J 3 . (3.17) 

give rise to the commutative subalgebra of the Sklyanin brackets (|1.3|) . Integrals of motion 
Ji,/ 2 and J3 are second, fourth and sixth order polynomials in variables Sj,tj. 

In variables ()2.7|) the Hamilton function 7T ()3.18|) is equal to I\ up to constants 
H = Ji + (a,b)(x- 2 C^-2(5 2 + <5 2 2 )) -2a0 = 

= ( J, AJ) + x" 1 (a x b, y x J) + 2(oa + ^b, J) (3.18) 

- x _1 (<5i - 5 2 )(a x b, y) - (<5i + <J 2 )(a x b, J), 
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where a = (ai, a 2 , £13), b = (b\, b 2 , 63) are numerical vectors and 

A = a®b + b®a, A^ = a^bj + a.,-^. 

Additional integral of motion 

Cl 



(3.19) 



K = 2x^ 2 - (a, b)(C 2 + 



4x 2 



(3.20) 



is a third order polynomial in momenta J. For brevity we present K at 8\ = 8 2 = only 



A" 



x 2 |J| 2 



x x (a x b, y x J) + 2(aa + f3b, J) 



(3.21) 



+ I J I 



8x 2 a/3 - 2(a, b) |y | 2 + 4x(aa - /3b, y) + (y x J, Ay x J 



- 4x(y, J)(aa — /3b, J). 
The third coefficient ^3 is a constant 



4x 2 2x 



4x 2 2x 



(3.22) 



Integrals of motion ()3.18|) and (|3.21J) depend on ten numerical parameters a*, 6j, a, /3, <5i, ^2 and 
they are defined up to canonical transformations. 

Suppose that Hamiltonian function has to depend on the third component of the vector 
y x J only. The Hamiltonian ()3.18j) has such form after rotation J — » UJ and y — > Uy with 
orthogonal matrix U defined by 



ei . e 2 




2 ' 








b = U- l b 



ei 




2 ' 



— 1 



e-2 




.0 



(3.23) 



where e« are eigenvalues of the matrix A ()3.19j) . Substituting a, b instead vectors a, b into 
()3.18|) and (|3.21j) one gets integrals of motion after rotation. 

If a and b are linearly dependent vectors, then a x b = and matrix A (J3.19)) has only one 
non-zero value e\ ~ |a| 2 . In this case the Hamilton function ()3.18|) 

E = 3\ + cJ x , ceR 

determines degenerate or superintegrable system with a noncommutative family of additional 
integrals of motion. For instance, there is the following quadratic integral 

A = ci J 2 + J 2 + Jf + c 2 y\ + c 3 (y 2 + y|) + c A y x J x + c 5 (y 2 J 3 - Vsh)- 

It is special case of the Poincare model 8J. For more details see jl()j . 
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If a x b 7^ 0, the matrix A (13.19)1 has two non-zero eigenvalues 

e i,2 = ( a i b) ± |a||b|, e 3 = 0. 
In this case after rotation and renormalization the Hamiltonian H ()3.18|) is equal to 

H 



H 



cJf-c 1 JI + k 1 (y 1 J 2 -y2Ji 



+ aJi + (3J 2 - (5 l + S 2 )J 3 - x \S 1 - 5 2 )y, 



(3.24) 



where 



c = — iy eie 2 1 , a — c\ — —iy 2e 2 1 (a + (3) 



P 



The Hamilton function H ()3.24)1 depends on five parameters instead of ten parameters in the 
initial Hamiltonian ()3.18jl . It allows us to impose constraints on the vectors a and b and to use 
triangular boundary matrices /C± [6[ instead of the general ones ()3.15J1 . The similar facts hold 
true for the BC n Toda lattices and for the Kowalevski-Goryachev-Chaplygin gyrostat jTT], 
which are also related with the reflection equations. 

Thus we can consider low triangular solution of the reflection equations /C + (|3.15|) 



b x - ib 2 







(3.25) 



without loss of generality. According to the Sklyanin method [3] the separated variables may 
be defined by entries of the following intermediate matrix 



T(.\): :/'(A)/C_(A) ( _° x J )7"(-A) 



1 
-1 



(3.26) 



which satisfies to the reflection equation [Hj 



r(\-v),T(\)T(v) 



+ T(A)r(A + v)T{v) -T{v)r{\ + v)T{\) . (3.27) 



(3.28) 



The separated coordinates q\ y2 are non-trivial zeros of the polynomial 

T 12 (A) = A(A 2 -^)(A 2 - ?2 2 ) = 0, 
whereas the conjugated momenta are equal to 

p k = -ikiT n (q k ) - ln(b 3 q k + a). (3.29) 

and the reflection equation 



We can prove that qk,Pk are Darboux variables using 
(|3~27l) . 
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Such us fC + is a low triangular matrix the generating function of integrals is 
r(A) = trace L(A) = -2(a, b)A 6 - h\ A - I 2 X 2 - h = 



= trace K+T(X) = (b 3 X + a)T u (\) + (h + ib 2 )\T 12 (\) - (b 3 X - a)T 22 (A) 
Substituting A = q k into this equation one gets two separated equations 

2(a, b) ql + hq k + h ql + h = exp(zp fc ) + det L(q k ) exp(-ip k ), k = 1,2. 

Here we took into account that Ti 2 (qk) = and T 22 (q k ) = det T(q k ) T{[ l {q k ). 

According to JT] we can introduce another separated variables related with the proposed 
separated variables by canonical transformation and by flip of parameters. Existence of the 
different separated variables is associated with the invariance of the Sklyanin brackets with 
respect to a matrix transposition T —* T*. 

4 Conclusion 

The Hamiltonians (|2.12j) and (|H.18j) belongs to the following class of the Hamiltonians 

H = (J, A3) + (a, y x J) + (b, J) + (c, y) , (4.30) 

possessing additional integrals of third and fourth degree in momenta. V.V. Sokolov kindly 
informed us that there exist three different such integrable cases only |12j . 

The research was partially supported by RFBR grant 02-01-00888. 

References 

[1] A.P. Veselov, Dokl. Akad. Nauk SSSR, v.270, p. 1094, 1983. 
A.I. Bobenko, Funct. Anal. Appl., v.20, p.53, 1986. 

[2] E.K. Sklyanin and T. Takebe, Comm.Math.Phys., v.204, p.17, 1999. 

[3] E.K. Sklyanin, J.Soviet.Math., v.31, p.3417, 1985. 

[4] V.B. Kuznetsov, A.V. Tsiganov, J.Phys.A., v.22, p.L73, 1989. 

[5] E.K. Sklyanin, Progr. Theor. Phys. SuppL, v. 118, p. 35, 1995. 

[6] I.V. Cherednik, Theor. Math. Phys., v.61, p.35, 1984. 
E.K. Sklyanin, J.Phys.A., v.21, p.2375, 1988. 

[7] V.V. Sokolov and A.V. Tsiganov, Theor. Math. Phys., v.131, p.118, 2002. 



8 



[8] H. Poincare, Bull. Astr., v.27, p.321, 1910. 

[9] H. J. de Vega and A Gonzales-Ruiz, J. Phys. A, v27, p. 6129, 1994. 
[10] A.V. Tsiganov, Reg. and Chaot. Dyn., v. 7, n.3, p.331, 2002. 
[11] A.V. Tsiganov, J. Phys. A., v.35, p.L309, 2002. 

[12] V.V. Sokolov, On a class of quadratic Hamiltonians on so(4), accepted to Dokl. Acad. 
Nauk, 2003. 



9 



